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Abstract 

We show that the use of wavelet bases for solving the momentum-space 
scattering integral equation leads to sparse matrices which can simplify 
the solution. Wavelet bases are applied to calculate the K-matrix for 
nucleon-nucleon scattering with the s-wave Malfliet-Tjon V potential. We 
introduce a new method, which uses special properties of the wavelets, for 
evaluating the singular part of the integral. Analysis of this test problem 
indicates that a significant reduction in computational size can be achieved 
for realistic few-body scattering problems. 


1 Introduction 

In this paper we show that bases of compactly supported wavelets can lead 
to a significant reduction in the size of the matrices used to solve few-body 
momentum-space scattering integral equations. When the kernel of the integral 
equation is expanded in a wavelet basis, the resulting matrix can be expressed 
as the sum of a sparse matrix and a matrix with small norm. Ignoring the small 
matrix leads to an accurate approximate solution. 

Our interest in wavelet bases is motivated by the successful application of 
wavelets to data compression problems. The coefficients of a wavelet basis pro¬ 
vide a compact representation for storing information that has smooth struc¬ 
tures on different scales. Wavelets are used for data compression in the FBI 
fingerprint archive and JPEG files |^. The property that makes wavelets 
useful for data compression suggests that they may also lead to sparse ma¬ 
trix representations of the kernels of the momentum-space scattering integral 
equations. Conventional methods lead to dense matrix representations of these 
equations; the computation of the solution of these equations for realistic models 
requires solving large systems of linear equations. 

While the use of wavelet bases lead to sparse matrices in both momentum 
space and configuration space, the resulting reduction in the size of the equations 
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in momentum space is significant, because the momentum-space kernels are 
represented by full matrices. These advantages are important for relativistic 
scattering problems, which are naturally formulated in momentum space. 

There are many different types of wavelets that are used in specific appli¬ 
cations. In this paper we use the compactly supported wavelets introduced 
by Daubechies [QQ in 1988 to evaluate the K-matrix for two spinless nucle¬ 
ons interacting with an s-wave Malfliet-Tjon V potential Q. Nucleon-nucleon 
scattering with the Malfliet-Tjon V potential provides a good test of numerical 
methods because the on-shell energy and strength of the attractive and repulsive 
parts of the potential lead to a problem with three different momentum scales. 

The Daubechies wavelets were the first non-trivial set of orthonormal wavelets 
of compact support discovered, and they have a number of properties that make 
them useful for numerical applications. However, the Daubechies wavelet basis 
functions involve a significant amount of structure that can make it difficult 
to appreciate some of these advantages. Some of the properties that we have 
found useful are briefly discussed below. They will be explained in more detail 
in section two. 

In applications, there are two bases, called the wavelet basis and the scaling 
basis. They are related by an orthogonal transformation, called the wavelet 
transform. This orthogonal transformation can be implemented efficiently using 
wavelet filters . The scaling basis is constructed from a single function, which 
is called the scaling function or father function. The scaling basis consists of 
discrete translates of a dyadic (power of 2) scale transformed copy of the original 
scaling function. The scale, which is called the fine scale, is determined by the 
application. The wavelet basis is also constructed from a single function, called 
the mother wavelet. The wavelet basis consists of discrete translations of dyadic 
scale-transformed copies of the mother wavelet function. These basis functions 
are called wavelets. This basis includes wavelets on all dyadic scales between 
the fine scale and a coarse scale, which is also dictated by the application. The 
wavelet basis is completed by translates of the scaling function on the coarse 
scale. 

The expansion coefficients for a function in the scaling basis are equiva¬ 
lent to weighted averages of the function over the support of the scaling basis 
functions. The wavelet basis provides an efficient representation to store this 
information on multiple scales. The principal advantages of the wavelet method 
are best explained in terms of properties of these two complementary sets of 
basis functions. 

1. The wavelets and scaling basis functions are orthonormal functions with 
compact support. This reduces global errors when these bases are used to 
represent local structures. 

2. The wavelet and scaling basis functions are related to a single “mother 
(resp. father) function” by unit translations and dyadic scale changes. 
This simplifies the computation of the basis functions. 

3. Finite linear combinations of the scaling basis functions can locally rep- 
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resent low-degree polynomials. This means that kernels which can be 
locally approximated by low-degree polynomials on the fine scale can be 
accurately represented by expansions in a scaling basis. This is the defin¬ 
ing property of the Daubechies wavelets. 

4. Wavelet basis functions are orthogonal to low-degree polynomials. This 
is the reason why the kernel in a wavelet basis is sparse. If the kernel 
can be locally approximated by a polynomial on a sufficiently small re¬ 
gion, the matrix elements corresponding to the wavelets with support on 
that region will be small or zero. In this representation the kernel is ac¬ 
curately represented by a small number of wavelet basis functions on a 
large scale. Since the wavelet basis is equivalent to the scaling basis, it 
follows that the wavelet basis can accurately represent kernels that can 
be locally approximated by low-degree polynomials on the fine scale by 
sparse matrices. 

5. All moments of wavelets or scaling basis functions can be computed exactly 
and efficiently. These moments are essential for performing the integrals 
needed for a typical calculation. 

6. There is a one-point quadrature rule which integrates the product of sec¬ 
ond degree local polynomials and scaling basis functions exactly. In our 
applications, the one-point quadrature rule reduces the Galerkin method, 
which constructs an approximate matrix equation by projecting the orig¬ 
inal equation on a subspace spanned by a set of basis functions, to the 
collocation method, which constructs an approximate matrix equation by 
demanding that the equation hold on a set of collocation points. The 
resulting systems of equations share the advantages of both methods. 

7. The wavelet transform provides a means to automatically determine the 
important basis functions. Specifically, it is enough to set all matrix ele¬ 
ments of the transformed kernel smaller than a specified minimal size to 
zero, which yields a sparse matrix. 

8. Integrals over the scattering cut can be computed accurately and efficiently 
using the scaling equation. These integrals can be expressed in terms 
of non-singular integrals that can be evaluated in terms of the known 
moments. 

9. The integral equation can be solved in the wavelet basis without computing 
a wavelet or a scaling basis function. This technique exploits many of the 
properties of the Daubechies wavelets and it saves the computation time 
involved in evaluating the basis functions. 

These features will be explained in more detail in the body of the paper. 

Given the many useful properties of Daubechies wavelets, it is natural to 
ask why they have not been extensively used to treat scattering problems. We 
speculate that part of the problem is the structure of the wavelets. Wavelets 
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are different than the smooth functions that are normally used in numerical 
analysis. They are finite linear combinations of solutions of a scaling equation, 
which relates a function on one scale to a linear combination of translates of the 
same function on a smaller scale. The solution of this scaling equation has a 
self-similar or fractal structure, which repeats itself on arbitrarily small scales. 
Functions with fractal structure are not amenable to numerical methods that 
exploit the smoothness of functions on a sufficiently small scale. Our experience 
is that typical adaptive quadrature methods are very inefficient when applied to 
wavelets. We found this to be especially true for integration over the scattering 
singularity. 

While the fractal structure of wavelets limits the utility of conventional nu¬ 
merical methods, the scaling equation provides an alternative tool, which when 
properly exploited, can be used to overcome all of the difficulties that arise with 
conventional methods. An important message in this paper is that success is 
obtained by embracing the scaling equation. 

In preparing this manuscript we found a number of references on wavelets 
to be useful in understanding how to apply wavelet bases to solve integral equa¬ 
tions. These include |H| H 111 H H FI] and Q. In addition, more details 

can be found in our web notes |13|, which provide a detailed summary of many 
of the things that we learned from these references. A key element that we 
believe is new in this paper is the treatment of the singular integral, which has 
proven to be both stable and accurate. This is also discussed in detail in these 
notes. 


2 Wavelets 

In this section we define Daubechies’ initial family of compactly supported or¬ 
thogonal functions Q. The defining property of Daubechies’ wavelets is that 
they are orthogonal to low-degree polynomials. 

As mentioned earlier, the basis functions consist of discrete translations and 
scale transforms of two functions, called the scaling function, which we denote by 
4>{x), and mother wavelet function, which we denote by ipix). An approximation 
space has two bases related by a fast orthogonal transformation, called the 
wavelet transform. One basis consists of translates of the scaling function on a 
fine scale. The second basis consists of translates of the scaling function on a 
coarse scale and the wavelets on all dyadic scales between the fine and coarse 
scale. The first representation will lead to a dense matrix representation of the 
integral equations while the second results in a sparse matrix representation. 
The choice of fine scale defines the approximation space; it is determined by the 
finest structure in the driving term and the kernel of the integral equation. 

We begin by discussing the scaling function, (j){x), which has compact sup¬ 
port and integrates to a value of one. To construct the scaling function we 
utilize two operations, dilation and translation. They are implemented by the 
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unitary operators: 


( 1 ) 


which is a dyadic scale transformation and 

f(l)(x) =(l)(x-l). (2) 

which is a unit translation. The scaling function (ji^x) is the solution of the 
scaling equation 

2iC-l 

Dcl^ix) = ^ (3) 

;=o 

with normalization 

/0(x)* = l. (4) 

Solutions, 4>{x), of the scaling equation with finite K are of special interest 
because these solutions have compact support in the interval [0, 2K — 1] [ p^ . 

The scaling equation is a linear renormalization group equation that re¬ 
lates the scaling function on a coarse scale to linear combinations of the scaling 
function on a fine scale. The coefficients hi, called the scaling coefficients, are 
numerical coefficients which define the type of scaling function. For most com¬ 
putations, a knowledge of the hz’s is all that is needed. In some wavelet literature 
the normalization convention for the scaling coefficients, hi, differs by a factor 
of V2 

The scaling basis functions, are defined by taking j dilations and k 

translations of the scaling function: 

</>,,fe(a:) := - k) ■ (5) 

Our convention, which follows refs. |4|Q, is that —oo < j < oo, with de¬ 
creasing j corresponding to finer scales and increasing j corresponds to coarser 
scales; some wavelet references use the opposite convention, with increasing j 
corresponding to hner scales. The scaling function cj){x) is defined so that the 
translates of the (j)j,kix) on a hxed scale j are orthonormal; 


= ^kl- ( 6 ) 

Self-consistency of the scaling equation (|^) and the orthonormality requirement 
(i constrain the coefficients hi so that 

2K-1 2K-1 

hi = V2, hihi-2n = SnO- (7) 

1=0 1=0 

For K = 1 these conditions fix the two coefficient /iq and hi uniquely. The 
K = 1 solution is called the Haar in scaling function. For larger values of 
K, additional constraints are needed to determine the coefficients hi. These 
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constraints will be discussed later in the paper. In all cases the, the number 
of hi ’s must be even and the support of the scaling function can be shown 
to be in the interval [0, 2K — 1]. The support of has width 2^{2K — 1). 

That the scaling equation has non-trivial solutions with compact support was 
not known before Daubechies seminal paper Q. 

The subspace of L^(R) spanned by 4’j,k{x) for fixed j is denoted by Vj. This 
space is the approximation space in our calculations. The parameter j defines 
the scale of these functions, and the size of the support of (f)j^k{x). The scale 
gets finer as j decreases, and j is determined by the finest scale of the problem. 
Two important properties of the different subspaces Vj are that 

V,_i D V, (8) 

and that Vj becomes dense in L^(R) as j —> — oo. Equation (||) is a consequence 
of the scaling equation. 

The mother wavelet "tpix) is defined to be orthogonal to the scaling function 
and its translates. The mother wavelet is associated with changes in the rep¬ 
resented function. It can be visualized as a small wave that integrates to zero. 
To define the mother wavelet recall Eq. (||). For each j define Wj to be the 
orthogonal complement of Vj in Vj-i: 

V,_1=W,©V,-. (9) 

Orthonormal basis functions of Wj are linear combinations of the scaling basis 
functions of V^-i. The different basis functions of Wj are discrete translations 
of a single function. The mother wavelet ip(x) is the element of Wo defined by 
the following linear combination of the scaling basis functions on V_i: 

2K-1 

ip(x) := b~^ ^ gkb4>{x) (10) 

1=0 

where the coefficients 

gi ■= (—(11) 

are fixed by the orthonormality conditions: 

(f V, T™V’) = 5im, (TV, T™V = 0. (12) 

The wavelets iIjo,i{x) ■= bipix) span Wo, while 

WV:=T^TV(V- (13) 

is an orthonormal basis for Wj. Note that the number of basis functions doubles 
as j is decreased by one. The support of ipj.iix) is identical to the support of 
4)jj(x)^ which is 

supp(V.i(a:)) = supp(V,i(a:)) = [2H,2\l + 2K - 1)]. 
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For a large negative value of j we can iterate to decompose Vj on a fine scale 
into a direct sum of Wfe’s on coarser scales: 

Vj = yVj+l © y^j+2 © • ■ ■ © Wj+m © Vj+m- (14) 

In wavelet terminology this is called the multi-resolution decomposition of Vj. 
The right- and left-hand sides of Eq. (Q) define two bases for the space Vj. The 
basis on the left consists of just the scaling basis functions (l)jj{x) on the fine scale 
“j”. The basis on the right consists of wavelets, ijjkjix), on all scales k from the 
finest (j + 1) to the coarsest (j + m) and the scaling basis functions 4>jj.rn,i(x) 
on the coarsest scale, (j + m). Note that unlike the spaces Vj, the spaces 
Wj on different scales are orthogonal. This leads to the following alternative 
representations of functions in Vj: 


j+m 

f{x) = ^ Cl(j)j,l{x) = ^ dl4>j+rn,l{x) + X! X! ^k,l'4’k,l(,x) (15) 

I I k—j-\-l I 

where the expansion coefficients are 

Cl ■■= {(l}j,i, f), di := {(j)j+jn,i,f) dk,i ■■= {ipk,i, f)- (16) 

It follows from Eqs. ( p^ and dil) that there is an orthogonal transformation 
relating the coefficients ci to di and dk,i- 

When we approximate a function f{x) as a linear combination the scaling 
basis functions, 4>jj{x), and the original function f{x) is approximately constant, 
f{x) = fi, on the support of (l)jj{x), the expansion coefficients, c; = {f,4’j,i) ~ 
are up to a scale factor, the value of the function at any point in the 
support of 4>jj{x). For a continuous function f{x), and a sufficiently fine scale 
j, it follows that, up to a scale factor, the coefficients ci track the value of / in 
the sense: 

Cl « 2^/V(2^0- (17) 

The coefficients c; form a very inefficient approximate representation of the 
function /(x); while the coefficients di and dkj form an efficient representation, 
since many of the dk,i will be small and can be neglected without introducing 
substantial errors. 

The Daubechies wavelets of order K are defined by the conditions that the 
mother wavelet satisfy 


J x^ip{x)dx = 0, 0<fc<iF — 1. (18) 

Using Eq. (^^ gives the additional constraints needed to determine the scaling 
coefficients hi: 


2K-1 2K-1 

0 = E 9^^'^ = E 0<k<K-l. (19) 

z=o z=o 
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Table 1: Scaling Coefficients 


hi 

K=1 

K=2 

K=3 

ho 

1/^2 

(1 + V3)/4y2 

(1 -b yiO -b \/5 -b 2yT0)/16V2 

hi 

1/V2 

{3 + i/3)/Ai/2 

(5 -b VTO + 3\/5 -b 2V^)/lQo/2 

h2 

0 

(3-v^)/4V2 

(10 - 2v^ -b 2^5 -b 2yiI))/16V2 

ho 

0 

(1 - v^)/4v^ 

(10 - 2v^ - 2\/5 -b 2yT0)/16V2 

hi 

0 

0 

(5 -b yiO - 3\/5 -b 2yi0)/16y2 

ho 

0 

0 

(l-b^lO- v'5-b2VT0)/16V2 


These equations, along with Eq. (^, determine the Daubechies scaling coeffi¬ 
cients, hi, of order if up to a reversal in the order of the coefficients. Note that 
the second equation in (Q) is actually several equations. There is a different 
equation for each value of n. For the case K = 1,2,3 these equations can be 
solved analytically. The solutions are given in Table 1 ^ . For larger values of K 
these equations can be solved numerically and tables of the scaling coefficients 
hi for K > 3 can be found in the literature §@- ft is straightforward to show 
that Eq. (^) holds on all scales: 


/ 


ipj^i{x)x^dx = 0 , — 1 . 


( 20 ) 


Plots of the K = 2 and K = 3 scaling function cj){x) and mother wavelet V'(^) 
are shown in Figures 1 and 2. The functions with K = 2 have support on the 
interval [0,3] while the K = 3 functions have support on [0,5]. The scaling 
function and mother wavelet can be distinguished by the conditions that the 
integral of the scaling function is one while the integral of the mother wavelet is 
zero. In both figures the compact support and fractal nature of these functions 
are apparent. 

One of the key features of the Daubechies wavelets follows from Eq. (p0|). 
Note that in the limit j —> —oo the expansion (IT^ becomes exact. If we expand 


3 

= '^di(j)j^i{x) + ^ y^^dk,ii’k,i{x), 

I k— — oo I 

then the condition (^) gives dk,i = 0 ior k < K, or 

i 


( 21 ) 


( 22 ) 


While a;^ is not a square-integrable function, at any point x the support condi¬ 
tions imply that the above sum has no more than 2K — 1 non-vanishing terms. 
It follows that on any compact subset this sum is pointwise identical to x^. Note 
that this is true for any scale j. This explains why the scaling function basis 
and the equivalent wavelet basis can be used to locally represent polynomials. 
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Eq. d^ ) explains why the second basis in Eq. is so efficient. If the 

function f{x) can be locally approximated by a polynomial of degree K on the 
support of tpj,iix), then the coefficients dk,i '■= i'4’k,i,f) ~ 0. 

Although the wavelets with higher values of K have some nicer properties, 
their support and the number of non-zero basis functions at any given point x 
increases with K. Thus, there is a point where the additional smoothness 
of a higher-order wavelet is balanced against the cost of manipulating functions 
with larger supports and more scaling coefficients, hi. For our application, we 
believe that the Daubechies K = 3 wavelet is the optimal choice to balance 
accuracy with computational speed. 

The important observation about the Daubechies wavelets is that once the 
scaling coefficients hi are determined, all of the calculations can be reduced to 
manipulations involving just these coefficients. This is done by making judicious 
use of the scaling equation to derive linear relationships for the quantities that 
we desire. This method is used to calculate all moments of the scaling and 
wavelet functions exactly and to construct quadrature rules. It can also be used 
to compute both tl’ix) and 4>{x) exactly at all dyadic points. Furthermore, the 
scaling equation leads to the transformation that we need to transform between 
the scaling basis function representation and the wavelet basis representation. 
The scaling equation is also used to accurately evaluate the integral of the 
product of the scaling basis functions and the scattering singularity. Specific 
applications of the scaling equation that are used to solve for the A'-matrix are 
outlined below. 

A. Computation of (j){x) and ^(x): 

To calculate </>(a;) we use the scaling equation in the form 

2K-1 

cj){x)= V2hi(l){2x - 1). (23) 

1=0 

Setting X to integer values n gives linear relations among the non-zero ^(n)’s at 
integer points: 

2K-1 

= E V2hi(j){2n-l). (24) 

1=0 

These homogeneous equations can be supplemented with the inhomogeneous 
equation 

n 

to get a linear system of equations for (p{x) at all integer points. Eq. (|5|) follows 
by using the normalization condition (^ to compute the expansion coefficients 
(© for f{x) = 1, and evaluating the result at an integer value of x. 

Given the values recursive use of the scaling equation ( p^ or (0) gives 
exact values of ipix) or ijj{x) at all dyadic rationale. Since these functions are 
continuous for A > 1, and the dyadic rationale are dense in the real numbers, 
(j){x) and tpix) can be calculated at all points by continuity. While this method 
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is exact at dyadic points, it is computationally intensive; however, the method 
discussed in this paper does not require the evaluation of the wavelets or the 
scaling basis functions. 

B. Moments: 

Moments of the scaling function are defined by 

< >ji := j (j)j,i{x)x^dx. (26) 

Using 

x^) = (</), f-^D-^x’^) (27) 

these moments can be expressed in terms of the moments < x^ >oo- Using the 
unitarity of the dilation operator and the scaling equation in 

< >00= (0, x'^) = {bcf), bxb (28) 


gives 


< x" >00= 


E 


fc! 


jk—m 


2^ ^ \/2 ^ ml(k — m)l 

I '' m—0 


< x"^ >00 


(29) 


This equation can be used to express the higher moments in terms of the lower 
moments 


2’^ - 1 


, k — 1 

e4e 


< a; >00= 

fc! 


-v/2 •>, m\{k — m)l 

'' m—0 ^ ^ 


< x'" >00 


(30) 


which can be used with the normalization condition (^ to recursively compute 
moments of any order. Given the moments of the scaling function, Eq. & can 
be used to compute moments of the wavelets. 

C. Quadrature Rules: 

Given a set of moments it is always possible to use them to construct quadra¬ 
ture rules. In most elementary applications a set of quadrature points {a;^} is 
given. The system of linear equations 


'^{xkTwk =< a:” >00 (31) 

k 


is used to solve for the weights Wk in terms of the < a;" >oo moments. For 
applications involving a scaling basis function on a fine scale only a small number 
of quadrature points are needed. 

In this paper most of the integrals are performed using a simple and effective 
one-point quadrature rule that follows from the identity 


< >00 = < >00 = 


2K-1 


71^"“ 


Z =1 


(32) 
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which holds for the moments of the Daubechies K > \ scaling basis functions. 
In this case one obtains 


(j){x)p{x)dx = p{< x^ >oo) for p{x) = a + bx + cx^ 


(33) 


which provides a one-point quadrature rulejll , with point x =< x^ >oo and 
weight w = 1, that integrates the product of the polynomial and the scaling 
function exactly. Translates and dyadic scale transformations of these moments 
give one-point quadrature rules for each of the 4>j^i{x). In applications we use 
points < x^ >ji and weights Wji for 4>j^i{x) on a small scale j: 


< x^ >ji= 2^< x^ >00 +0 Wji = 2^1"^ 


(34) 


This quadrature rule is used on the finest scale with the Daubechies K=2,3 
scaling basis functions, which can represent local polynomials of degree one and 
two exactly. For Daubechies wavelets with larger K-values, which can locally 
represent polynomials of degree iF — 1, multi-point quadratures rules can be 
constructed which exactly integrate the product of the scaling function with 
degree K — \ polynomials. 

D. Partial Moments: 

In solving integral equations on a finite or semi-infinite interval, it is neces¬ 
sary to calculate integrals where an end-point of the interval is in the support 
of a scaling basis function § 0 - By using the scaling and dilation operators 
these integrals are simply related to the following integrals: 

poo /*0 

:= / f'^(j){x)x'^dx, I^_ := / f'^(j){x)x'^dx (35) 

and 

poo pO 

Cn+ ■■= / f"^H^)f-mdx, := / f-^^ix)f-^ix)dx. (36) 

Jo J — oo 

To calculate these quantities one can use the scaling equation in the form (^. 
This leads to the linear relations 


rk _ cx — k—Yjl 

^m-\- ^ 


2K-1 

^ hil!: 




Z=0 


(37) 


2K-1 2K-1 

^mn-\- — E E hj-hsl2m-\-r,2n-\-s,-\- (38) 

r—0 s—0 

Because I^_^_ =< x^ >om for m > 0, Eq. (|3|) becomes a linear system for the 
unknown partial moments in terms of the full moments. These equations can 
be solved for the non-trivial 7^+. The are obtained using 
Useful examples can be found in [|^ and Q . 
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For Imn-\- we have for m > 0 or n > 0: 

^mn-\- ^mn- 


( 39 ) 


from the orthonormality of the (j)oi (cc). This leads to a small linear system that 
relates the unknown integrals Imn+ to the known Imn = 5mn- The Imn- can 
again be obtained by subtraction. 

These linear systems can be solved for the non-trivial Imn+ ■ The correspond¬ 
ing partial moments on other scales can be obtained from the Imn+ using the 
unitarity of the operator D on the half intervals. Partial moments with different 
endpoints can be obtained by translation and subtraction. These integrals are 
used to treat endpoint quadratures after the if-matrix equation is transformed 
to a finite interval. 

E. Singularity: 

The scaling equation can also be used to solve for the singular integrals of 
the form 



(j){x — k) 
X ± i0+ 


dx. 


(40) 


Application of the scaling equation (||) leads to the linear relations: 


2K-1 


1=0 


(41) 


Using this equation for different values of k leads to a set of linear equations that 
relate the for values of k where the scaling basis function has support con¬ 
taining the singularity to the for values of k where the scaling basis function 
has support far away from the singularity. One more equation relating theses 
quantities is needed to pick up the contribution from singularity. Specifically 
we use 

n 

which follows from Eq. (^^ to get: 


ITT = 


dx 


X ± i0+ 




dx 


X ± i0+ 


= 

/ ^ -^n:a ’ 


(43) 


where 



(p{x — n) 

X ± i0+ 


dx = 


(44) 


except when the points ±a are in the support of (f>nix). The values a is chosen 
to be far from the singularity. For the principal value integral, which is used in 
this paper, we replace ^iir on the left of Eq. ( ^ by zero. 

The input needed to solve this system are the values of and far from 
the singularity. For |n| sufficiently large we can calculate these quantities by 
expanding them in terms of the moments 



>00 


(45) 
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Table 3 - Integrals over singularity: K.—2 



-2.779949550280 

Tz 4.291495373146 

1-2 

-0.269952669589 

±i 1.149902719556 


Table 4 - Integrals over singularity: K=3 



-0.1717835441734 

=Fz 4.041140804162 


-1.7516314066967 

±i 1.212142562305 

1-3 

-0.3025942645356 

=Fz 0.299291822651 

I-A 

-0.3076858066180 

Tz 0.013302589081 


and the partial moments 



— n) 1 
x±iQ+ ^ ^ n 




< X 


>00:a ■ 


(46) 

These integrals can be accurately approximated for sufficiently large |n| by trun¬ 
cating the series, ( ^ ) and (^, for sufficiently large |n|. 

Examples of the integrals that overlap the singularity are given in Tables 3 
and 4 for the K = 2 and K = 3 Daubechies wavelets. The endpoint integrals, 
where the singularity is located at the boundary of the support, are not singular 
because of the continuity of (j){x) for K > 1. These integrals, and integrals over 
scaling basis functions close to the singularity are computed using the linear 
equations © and Singular integrals corresponding to different scales are 
related by the identity 


/ 


(j)j^k{x)dx 

X ± j0+ 




(47) 


The principal-value integrals are given by the real parts of 


3 Test Problem 


The power of the wavelet method is demonstrated by calculating the 7f-matrix 
for an s-wave Malfliet-Tjon V potential. The momentum-space integral equation 
for the two-body 7f-matrix is 

Ko{pi,P2,Po) = vo{pi,P2) - mP.V. f ’^^2 ^oip',P2,Po)dp' (48) 

Jo P — Po 

where m is the nucleon mass and P.V. indicates the principal value. The inter¬ 
action potential is a sum of two Yukawa interactions 


r ^ 1 f Pi + Pi+P2 + 2piP2 ^ 

Vo{Pi,P2) ■■= y-z -In ———2 -^- 

^ 2tTPiP2 \pi +Pi+Pi- 2piP2 J 

with strength and range parameters given in Table 2. 


( 49 ) 
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Table 2: Potential Parameters 


Ijin 

Ai 

Ml 

A 2 

M2 

41.47 MeV fm^ 

-570.316 MeV fm 

1.55 fm-^ 

1438.4812 MeV fm 

3.11 fm“^ 


Test calculations are done for the half on-shell if-matrix, Kq(pi,P 2 ,po), with 
P 2 =Po, and pl/m = 10 and 80 MeV. 


4 Wavelet Techniques for Singular Integral Equa¬ 
tions 


We solve the integral equation by transforming the half-interval [0,oo) to 
a finite interval [—a, b] with the singularity po transformed to the origin. 

For fixed values of p 2 and po define 


f{pi):=K{pi,p2,po) (50) 

g{pi) :=vo{pi,p 2 ). (51) 

The fixed variables p 2 andpoi which are just parameters in the integral equation, 
are suppressed in this notation. 

If we also define the non-singular part of the kernel 


L{pi,P2) 


voiPi,P2)pl 

m - 

P 2 +Po 


(52) 


the integral equation for the if-matrix has the form 


f{Pi) = gipi) - P-V. [ f{p')dp'. 

Jo P - Po 


(53) 


To transform [0,oo) to the interval [—a,b] we use the following map, which 
also maps the singularity at po to zero: 


Let 


, , b a-\-u 

P=P[u)\=pq u = u(p) 

a 0 — u 

_ ab{p-po) 
ap + pob 

(54) 

b {b +a) 

dp = P{j ,, .^du. 

a[b— uy 


(55) 

1 a{b — u) 1 

p- Po (ad- b)po u 


(56) 

. / b a + u\ 

g{u):=g[po 

\ a b — u J 

•> 

(57) 
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and 


/(w) := / 


b a + u\ 

Po-T - 

a b — u J 


(58) 


, , , b a + u b a + V 

L{u,v) = L po-- - ,Po-T - 

' a b — u a b — V 


Po 


(b-v)’ 


(59) 


where the factor b/{b— v) is the product of the Jacobian and the residue of the 
transformed singularity 


b — V a (b — vY (a + b)pQ 

The integral equation in the u, v variables has the form 


b (b + a) a{b — v) 


f{u) = g{u) + 


L{u, v) 


f{v)dv. 


(60) 


(61) 


To solve this we represent f{u) using the scaling basis on a sufficiently Hne scale 

j- 

/(m) « ^/„^j,„(w). (62) 


Inserting the expansion (|^ in the integral equation (^^ gives 


=giu)+Y^ 


L{u,v] 


4>pri.{u)dv fn- 


(63) 


To get a linear system for the coefficients /„ we multiply the above equation by 
and integrate from —a to b: 


^ ^ d^mnfn — g-m + ^ ^ 


rb L(u,v) 

- —(j)j^niu)dudvfn, (64) 


where ^ 

Nmn ■= / 4>pm{u)4>pn{u)du. (65) 

J —a 

The matrix Nmn is Jmn except when the support of (pj^rn{u) and 4>j^n{u) contains 
—aoib. This is expressed in the form 

Nmn — dmn T ^mn- (66) 

The integrals Nmn can be computed exactly using Eq. (|3^^. The integrals, 
Nmn, the singular integrals, as well as the moments, partial moments, and 
quadrature points can be calculated once and stored for later use. 

The vector gm is given by 


5m = / g{u)(j)j^m{u)du, (67) 

J —a 
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which is computed using the one-point quadrature when —a or b is not in the 
support of 4‘j,m(u). In the case that —a or 6 is in the support of we 
use the partial moments (|^) to construct endpoint quadrature rules (^). We 
use K 1 transformed Gauss-Legendre points. The weights are determined by 
solving a linear equation determined by the requirement that the quadrature 
rule reproduces the lowest-order partial moments up to order K. 

The matrix elements of the kernel are evaluated using 


T™, . — 


4^j,m 


L(v, u) 


du dv 


— / 


L{v, u) — L(v, 0) 


4>j^n {u) du dv 


+ [ (j)j^rn{v)L{v,0)dv ( du. (68) 

J — a J— a ^ 

For most basis functions the v integral is done using the one-point quadrature 
rule. This is also the case in Eq. (O) for the “u” integral and for the integral 
over the subtracted term in Eq. (|68[). The quadratures used for the endpoints 
of gm are used to integrate scaling functions with support on the endpoints of 
the integral. The singular “u” integral on the last line of Eq. is done using 
the method discussed in the previous section. 

The resulting equation has the form 


fm — 9 m ri 


^m,n)fn- 


(69) 


The coefficients fm obtained by solving this matrix equation can be used to cal¬ 
culate a refined solution, which is constructed by substituting the series solution 
into the right-hand side of the integral equation. If the integrals are evaluated 
following the same procedure used in the integral equations, the solution f{u) 
can be obtained for any value of u without having to calculate the basis func¬ 
tions 4’j,m{u). This is signihcant because the evaluation of the basis functions is 
computationally intensive. The refined solution, /(m), can be used to evaluate 
the half on-shell if-matrix at any point p by using the inverse of (^8|): 

K{p,Po,Po) = f{p) = f (70) 

V ap + pob ) 

The calculation outlined above, where the transformed K-matrix equation is 
solved in the scaling function basis, has a full kernel matrix. While the treatment 
of the singular integral is very stable and accurate, the method does not lead to 
a reduction in the size of the matrix equations when compared to calculations 
using splines or other numerical techniques. In the next section we show that if 
Eq. ( 1^ is transformed to an equivalent equation in the wavelet basis, that the 
resulting kernel matrix is well approximated by a sparse matrix, which leads to 
a considerable reduction in the scale of the numerical calculation. We also show 
that the wavelet transform can be implemented by an efficient algorithm that 
uses a simple filter on multiple scales. 
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5 Wavelet Transform 


The computation done in the previous section leads to a linear system with a 
large dense matrix. The solution is a set of numerical coefficients fi = Cj^i of 
the approximate solution in the scaling function basis for a fine scale j. As 
mentioned earlier, there are two equivalent representations for the solution /(u) 
given in Eq. (H- The wavelet transform for a problem with 2^ basis functions 
maps coefficients of the basis functions on model space Vj to coefficients of the 
basis functions on the equivalent space Wj+i © • • • Wj+N © Vj+N- 

The form of the transformation is determined by applying the Eqs. i) and 
to map 

Vi ^ Vi+i (71) 

and 

Vi ^ Wi+i. 

This is done N times, successively applying Eqs. and ( 0 ) toVj, Vj+i 
The coefficients on the next (coarser) scale are obtained using 


1 

(73) 


(74) 


i 


Because only a finite number of coefficients are considered, there is some ambi¬ 
guity in how to treat the coefficients that appear on the right-hand side of these 
equations, which are out of range. The simplest method to treat these terms is 
to use a periodic wrap-around condition: 


(72) 

■ Vj+N-l- 


^ Cj+k,i , (75) 

which is applied at each level, k. Since the result is an orthogonal transformation 
of the original basis, this choice does not affect the solution. It will have a small 
edge effect on the structure of the resulting sparse matrix. 

The structure of the successive transformations is illustrated for the case 
j = —3 and iV = 3 by the following diagram. In this diagram Cj^m represents 
the coefficients of the basis functions in Vj and dj^m represents the coefficients 
of the basis functions in Wj. The successive transformations act on C- 3^1 as 
follows: 


C- 3,1 \ 




( ^ 


( C-0,1 ^ 

C-3,2 


c-2,2 


C-1,2 


d-0,1 

C- 3.3 


c-2,3 


d-1,1 


d-1,1 

c-3,4 


c-2,4 


d-1,2 


d-1,2 

c-3,5 


d-2,1 


d-2,1 


d-2,1 

c-3,6 


d-2,2 


d-2,2 


d-2,2 

c-3,7 


d-2,3 


d-2,3 


d-2,3 

c-3,8 / 


\ d-2,i / 


\ d-2.4 / 


\ d-2,4 / 
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This transform can be implemented efficiently. Since it is an orthogonal trans¬ 
formation, if it is expressed in matrix form, the inverse is just the transpose. 
Both the wavelet transform and its inverse can be implemented in 0{N) |Q 
steps. 

The advantage of using the wavelet basis is that the coefficients dj^k in Eq. 
(15) will be small if / f{u)ipj^k{u)du « 0. This will be true, by virtue of Eq. 
(20), whenever f{u) can be locally approximated by a polynomial of degree less 
than K on the support of tpj^kiu). A practical rule of thumb states that the 
number of non-zero matrix elements grows like 10A^log(l/e) Q, where N is the 
size of one dimension of the matrix, and e is the absolute value of the ratio of 
size of the smallest matrix element retained to the largest matrix element. This 
behavior is consistent with our observations. 

In Figure 3 a simple example is used to illustrate how structure appears on 
different scales once the wavelet transform is applied. In this case a product 
of two Gaussian functions with displaced broad and narrow peaks is expressed 
in a K=3 Daubechies scaling basis. The result is transformed with the wavelet 
transform and the strength of the coefficients in the wavelet basis are studied. 
The figure shows the coefficients of the scaling functions on the coarse scale 
give a low-resolution large-scale representation of the original function. The 
coefficients of the wavelets on different scales builds the fine structure. Note the 
that scales on the vertical axes differ in each frame. It is also interesting to see 
how the wavelet transform automatically finds the parts of the original function 
with structure. 

To fully use the wavelet transform, the number of basis functions should be 
a power of two. In this paper the choice of —a and b is adjusted so that the 
number of basis functions is a power of two. By translating the interval [—a, 6] it 
is also possible to control the number of basis functions with support on either 
side of the singularity at zero. 

In applications the wavelet transformation is applied to the driving term 
and to the rows and columns of the kernel of the integral equation. This leads 
to an equivalent linear system with respect to the transformed basis. In this 
new basis we expect the matrix of the kernel to be the sum of a sparse matrix 
and a matrix with small norm. Ignoring the small normed matrix leads to an 
approximate linear system with a sparse kernel. 

This approximate solution can be transformed back to the scaling represen¬ 
tation and substituted back into the integral equation. Evaluating the integrals 
using the one-point quadrature and methods discussed previously, gives the 
refined solution, constructed without computing the wavelets or scaling basis 
functions. 


6 Results 

Initial calculations are performed using K=2 and K=3 Daubechies wavelets 
using N scaling basis functions, up to Nmax = 2® = 512. The finest resolution 
corresponds to j := J = —7. The values of a and b are given by a = 1 and 
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b = —a + {N — 2K + 2)2'’. The solutions are represented directly in terms of 
a series expansion in the scaling function basis and using this expansion in the 
integral equation to construct the refined solution. The refined solution has the 
advantages that the fractal structure of the basis does not appear in the solution 
and the basis functions do not have to be evaluated. These calculations were 
checked against conventional differential and integral equation programs, and 
the agreement was excellent. 

Calculations were done for the half on-shell K-matrix at 10 and 80 MeV. 
Tables 5-8 show the values of the refined on-shell K-matrix elements at 10 and 
80 MeV for the Daubechies K = 2 and K = 3 scaling function bases for different 
resolutions j (approximation spaces Vj). These tables illustrate the convergence 
of a full calculation as the resolution is increased. In each successive calculation 
the resolution is increased by a factor of 2. 

The convergence for the K = 2 and K = 3 calculations are compared for 
different values of j = J, which defines the smallest scale, in Tables 9-10. It 
is apparent that both methods are accurate, but the K = 3 Daubechies basis 
gives better results for the same number of basis functions. 

In this paper the refined solution with N = 2^ scaling basis functions serves 
as the benchmark calculation that is used to test the approximate wavelet cal¬ 
culations with sparse matrices. In general, the method leads to accurate and 
stable solutions. This solution is converged to the accuracy expressed in the 
tables and agrees with calculations using other methods. 

Applying the wavelet transform to the integral equation and throwing away 
all matrix elements of the kernel smaller than a given size leads to approximate 
equations with sparse matrices. Tables II-I3 compare the solution of different 
sparse matrix approximations to the full 512 basis function calculation for the 
K = 2 and K = 3 Daubechies basis functions at 10 and 80 MeV. The first 
column in each of these tables indicates that matrix elements with an absolute 
value smaller than e times the absolute value of the largest matrix elements are 
set to zero. The second column indicates the percentage of the 512 matrix ele¬ 
ments that are non-zero in the sparse matrix approximation. The third column 
indicates the on-shell value of the AT-matrix obtained by applying the integral 
equation refinement to the sparse matrix solution. Column five shows the rela¬ 
tive error at the on-shell point and column six shows the mean square error of 
solution to the integral equation over the whole domain defined by: 


/(o) - m 

/(O) 


(77) 


and 

lifju) - Mu))‘^du)’/^ 

(//2(u)du)i/2 ■ 

The power of the wavelet method is apparent in the K = 3 Daubechies basis. 
Tables 12 and 14 show, for the case e = 10“®, that about 96% of the matrix 
elements are eliminated relative to the full calculations. The on-shell error is a 
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few parts in 10^ for both energies. Even more impressive is the fact that the 
mean-square error for the half-shell if-matrix is also about the same order of 
magnitude. 

Figure 4 shows the transformed half on-shell K-matrix, f{u) at 80 MeV 
using the sparse matrix approximation with e = 10“®. This figure is a plot of 
the series solution, which is indistinguishable from the refined solution. The 
smoothness of the series solution is apparent from this figure. While in Figure 
5 we show the half on-shell K-matrix, f{pi), at 80 MeV using the sparse matrix 
approximation with e = 10“®. This plot shows the refined solution, which 
is obtained by substituting the series solution back in the integral equation. 
As in the transformed case, this is indistinguishable from the series solution. 
Figure 6 provides a graphic illustration of how the non-zero terms of the sparse 
kernel are distributed for the case e = 10“®. What is relevant is that the two 
dimensional plot is dominated by a number of one-dimensional structures. The 
vertical bands in this matrix correspond to basis functions with support near 
the singularity, where we expect significant structure. The different diagonal 
bands arise because local structures in the kernel can couple wavelets on several 
scales. 

While it is possible to calculate perturbative corrections to approximate 
solutions, it make more sense calculate the matrix equation in the scaling basis at 
a finer scale, and use the wavelet transform to get a better approximation. This 
is because the matrix with small norm, that is needed to compute perturbative 
corrections, is not sparse. 

The tables and calculations indicate that wavelet methods provide a poten¬ 
tially powerful technique for reducing the size of momentum-space scattering 
integral equations. In this paper we did this by transforming a dense kernel 
using the wavelet transform. These operations can be done very efficiently, re¬ 
sulting in a matrix which is the sum of a sparse matrix and a matrix with small 
norm. 

Solving the linear system using only the sparse matrix leads to accurate 
solutions. In our application a decrease in the size of the matrix by a factor of 
20 led to an error of a few parts in 10® for the Daubechies K = 3 wavelets. 

The speed of our computations is helped by using the one-point quadrature, 
which leads to accurate answers when used with the Daubechies AT = 2 or 3 
wavelets. Another nice feature of the wavelet method is that it automatically 
puts basis functions at places where the kernel has structure, as compared to 
splines, where the placement of the knots must be done by hand. 

7 Summary 

In this paper we adapted numerical methods based on wavelets to calculate the 
K-matrix for a Malfliet-Tjon V potential in momentum space. The motivation 
for this work was to determine if numerical methods based on wavelets could be 
used to reduce the size of the matrices that occur in numerical implementations 
of the momentum-space integral scattering equations. 
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There are three reasons for investigating this problem: 

1. Relativistic formulations of the few-body problem are naturally formu¬ 
lated in momentum space. This is because the momentum operators are 
inhnitesimal generators of the Poincare group. Operators, like Wigner 
rotations, have simple forms in momentum space. 

2. Boundary conditions and scattering asymptotic conditions can be cleanly 
formulated in momentum-space compact-kernel integral equations. 

3. The kernel matrix that appears in most numerical implementations of 
momentum-space scattering equations is full. The size of the kernel matrix 
is also large for few-body calculations with realistic interactions. 

We found that the matrix representation of the momentum-space scattering 
integral equations in the wavelet basis leads to a linear system with a sparse 
approximate kernel matrix. This leads to a signihcant reduction in the size of 
the dynamical calculation. 

Our applications used the Daubechies K = 2 oy K = 3 wavelets. For each 
value of K there were two bases, the scaling function basis and the wavelet 
basis, which were related by a fast orthogonal transform, called the wavelet 
transform. The K = 2 bases can locally represent lines of arbitrary slope on 
subintervals corresponding to the finest scale. The K = 3 bases can locally rep¬ 
resent quadratic curves on subintervals corresponding to the finest scale. Our 
observations indicate that we obtain better accuracy with the same number of 
basis functions using the K = 3 basis than with the K = 2 basis without a sig¬ 
nificant increase in computational complexity. While it is possible to use a basis 
with a higher K value, in order to take advantage of the additional capabilities 
of the basis functions, the one-point quadrature would have to be replaced by a 
multi-point quadrature. The additional benefit is not worth the complexity and 
increased computational effort required by the additional quadrature points. 

One potential difficulty with using wavelet methods in numerical calculations 
is the self-similar structure of the basis functions. Functions that do not become 
smooth on a sufficiently small scale are not amenable to conventional numerical 
methods. Initially, we found that this property caused numerical challenges. 
To overcome these problems we utilized the scaling equation to exactly calcu¬ 
late moments of the scaling basis functions or a wavelets times a polynomial. 
We found that by using these moments and the scaling equation we were able 
to exactly compute overlap integrals, to accurately evaluate integrals over the 
scattering singularity, and to exactly compute endpoint overlap integrals. It 
is only after using the scaling equation to do all of the required quadratures 
that we were able to obtain stable and accurate answers. While many of the 
techniques discussed in this paper are standard wavelet procedures, the method 
for treating the scattering singularity seems to be new to this paper. 

As a basis for solving integral equations, the scaling function basis has many 
of the same advantages as a spline basis. The most important common features 
are the compact support and the ability to exactly represent polynomials on 
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small scales. This means the they can be used to efficiently approximate local 
structures. The scaling function basis has the additional feature that the func¬ 
tions are orthonormal. In addition, for the K = 2,3 Daubechies wavelets, the 
one-point quadrature can be used to exactly compute the integral of the product 
of a degree two polynomial and the scaling basis function on the support of the 
scaling basis function. By employing the one-point quadrature rule, the matrix 
formulation of the integral equation can be computed very quickly. While the 
wavelets and scaling basis functions are computationally intensive to compute, 
the one-point quadrature can be used to replace the evaluations of these func¬ 
tions if the series solutions are substituted back in the integral equation. This 
method, leads to our refined solution. The only input that is needed are the 
scaling coefficients, hi, which are known analytically for the K = 2 and K = 3 
wavelets. 

In working on this problem, we were able to find a number of previous ap¬ 
plications of wavelet methods to integral equations; but were unable hnd appli¬ 
cations that treated integral equations with the singular kernels of momentum 
space scattering equations. We found that the required integrals could be accu¬ 
rately computed using the scaling equation. The method has proved to be very 
stable. It has the advantage that it reduced the calculation of these integrals 
to solving a small system of linear equations, the results of which can be stored 
for later calculations. 

The real advantage of the wavelet method over the conventional methods 
occurs after the wavelet transform is applied to the K-matrix equation in the 
scaling function basis. The kernel matrix in the wavelet basis can be expressed 
as the sum of sparse matrix and a matrix of small norm. Ignoring the small 
matrix leads to an approximate solution. The case featured in Figure 6, used 
Daubechies K = 3 wavelets. The sparse matrix approximation used in Figure 
6 led to a 96% reduction in the number of non-zero matrix elements. This was 
done with a dimensionless mean square error of a few parts in 10® 

The beauty of the wavelet transform is that it is fast, and it automatically 
finds the basis functions needed to represent the small scale structure of the 
solution. 

In this paper we have exploited the scaling equation to formulate a powerful 
method for solving momentum space scattering integral equations. While we 
limited our application to what was necessary to solve the scattering integral 
equations, it is apparent that there are many thing that can be computed exactly 
or to a very high accuracy with the scaling equation. For example, it is possible 
to exactly integrate products of several scaling functions and powers of the 
variable. This suggests that these methods could be very useful for solving a 
large class of physics problems. 

Our conclusion is that wavelet bases provide a powerful method for solving 
momentum-space scattering integral equations. The required matrix elements 
are easy to compute, and the wavelet transform can be implemented very effi¬ 
ciently, leading to an approximate linear system with a sparse matrix. 
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Table 5 - Dependence on number of basis functions: K=2, E= 10 
MeV 


-J 

N 

series on-shell 

refined on-shell 

3 

32 

-124.681226 

-124.401416 

4 

64 

-124.924473 

-124.853374 

5 

128 

-124.984907 

-124.967026 

6 

256 

-124.999853 

-124.995372 

7 

512 

-125.003567 

-125.002445 


Table 6 - Dependence on number of basis functions: K=3, E=10 MeV 


El 

N 

series on-shell 

refined on-shell 

3 

32 

-125.051451 

-125.034060 

4 

64 

-125.007967 

-125.006049 

5 

128 

-125.005171 

-125.004948 

6 

256 

-125.004847 

-125.004820 

7 

512 

-125.004806 

-125.004803 


Table 7 - Dependence on number of basis functions: K=2, E=80 MeV 


El 

N 

series on-shell 

refined on-shell 

3 

32 

-6.53375948 

-6.38393342 

4 

64 

-6.45483277 

-6.41711946 

5 

128 

-6.43490787 

-6.42555390 

6 

256 

-6.42998750 

-6.42766546 

7 

512 

-6.42877076 

-6.42819277 


Table 8 - Dependence on number of basis functions: K=3, E=80 MeV 


El 

N 

series on-shell 

refined on shell 

3 

32 

-6.44161445 

-6.43154124 

4 

64 

-6.42926712 

-6.42868443 

5 

128 

-6.42842366 

-6.42840177 

6 

256 

-6.42837147 

-6.42837210 

7 

512 

-6.42836848 

-6.42836877 


Table 9 - Refined on-shell values, E=10 MeV 


N 

K=2 on-shell 

K=3 on-shell 

32 

-124.401416 

-125.034060 

64 

-124.853374 

-125.006049 

128 

-124.967026 

-125.004948 

256 

-124.995372 

-125.004820 

512 

-125.002445 

-125.004803 
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Table 10 - Refined on-shell values, E=80 MeV 


N 

K=2 on-shell 

K=3 on-shell 

32 

-6.38393342 

-6.43154124 

64 

-6.41711946 

-6.42868443 

128 

-6.42555390 

-6.42840177 

256 

-6.42766546 

-6.42837210 

512 

-6.42819277 

-6.42836877 


Table 11 - Sparse Matrix Convergence: K.—2, E=10 MeV, 3—-7 


e 

percent 

on-shell value 

on-shell error 

mean-square error 

0 

100 

-125.00245 

0 

0 

10 “® 

34.24 

-125.00244 

6.94x10-® 

8.76x10-® 

10 “® 

20.8 

-125.00244 

3.82x10-® 

7.81x10-® 

10-7 

12.16 

-125.00440 

1.57x10-® 

1.63x10-® 

10 -® 

6 

-124.99271 

7.78x10-® 

9.53x10-® 

10 “® 

2.61 

-124.96069 

.000334 

.000435 

10-4 

1.16 

-124.47119 

.00425 

.00457 

10 -® 

.55 

-120.42759 

.0366 

.0415 

10-2 

.31 

-116.14345 

.0709 

.154 

Table 12 - Sparse Matrix Convergence: K= 

:3, E=10 MeV, J= 

e 

percent 

on-shell value 

on-shell error 

mean-square error 

0 

100 

-125.00480 

0 

0 

10 -® 

17.78 

-125.00480 

1.05x10-® 

2.56x10-® 

10 -® 

11.38 

-125.00480 

5.14x10-® 

2.44x10-7 

10-7 

6.6 

-125.00475 

4.49x10-7 

1 .88x10-® 

10 -® 

3.76 

-125.00269 

1.69x10-® 

2.08x10-® 

10 -® 

2.14 

-124.99030 

.000116 

.000228 

10-4 

1.24 

-124.85112 

.00123 

.00217 

10 -® 

.72 

-123.82508 

.00944 

.0117 

10-2 

.38 

-125.25766 

.00202 

.128 


Table 13 - Sparse Matrix Convergence: K.—2, E= 80, MeV 3—-7 


e 

percent 

on-shell value 

on-shell error 

mean-square error 

0 

100 

-6.4281928 

0 

0 

10 -® 

38.25 

-6.4281928 

4.16x10-17 

6.30x10-® 

10 -® 

23.59 

-6.4281927 

1.54x10-® 

1.13x10-7 

10-7 

13.84 

-6.4281936 

1.30x10-7 

2 .11x10-® 

10 -® 

6.89 

-6.4282684 

1.18x10-® 

3.87x10-® 

10 -® 

2.91 

-6.4286669 

7.38x10-® 

.000211 

10-4 

1.18 

-6.4286269 

6.75x10-® 

.00166 

10 -® 

.55 

-6.4083948 

.00308 

.0113 

10 -® 

.3 

-6.3206920 

.0167 

.101 
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Table 14 - Sparse Matrix Convergence: K=3, E=80, MeV 3—-7 


e 

percent 

on-shell value 

on-shell error 

mean-square error 

0 

100 

-6.4283688 

0 

0 

10 “® 

19.99 

-6.4283688 

1.52x10-1° 

1 .20x10-® 

10 “® 

12.94 

-6.4283690 

3.44x10-® 

2.06x10-1' 

10 -'^ 

7.42 

-6.4283703 

2.33x10-1' 

1.87x10-° 

10 “® 

4.08 

-6.4283333 

5.51x10-° 

4.38x10-° 

10 “® 

2.22 

-6.4278663 

7.82x10-° 

.000994 

10-4 

1.21 

-6.4244211 

.000614 

.00845 

10 “® 

.67 

-6.4154328 

.00201 

.0229 

10-2 

.34 

-6.2935398 

.021 

.102 
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Figure 1: Daubechies K = 2 Scaling function and mother wavelet 
function 
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Figure 2: Daubechies K = 3 Scaling function and mother wavelet 
function 
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Figure 4: Transformed K-Matrix - E = 80 MeV - expansion in wavelet 
basis 
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Figure 5: K-Matrix - E = 80 MeV - refined solution 
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Figure 6: Sparse Kernel for K = 3, e = 10 ® 





